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DISCLAIMEP 



Although each program has been tested by Its contributor, no 
warranty, express or Implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USEBS Group, in con- 
nection therewith. 
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Program No, _____________ Date 

Program Name: 



1. Does the abstract adequately describe what the program is and what Yes No 

it does? 



Comment 



2. Does the program do what the abstract says ? Yes No_ 

Commen t 

3. Is the Description clear, understandable, and adequate? Yes No 

Commen t 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No. 

Commen t 

Are the Sense Switch options adequately described (if applicable)? Ye s No 

Are the mnemonic labels identified or sufficiently understandable? Yes No" 
Commen t """""" 

5. Does the source program compile satisfactorily (if applicable)? Ye s No 

Commen t 

6. Does the object program run satisfactorily? Yes No 

Commen t 

7. Number of test cases run , Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

Commen t 

8. Does the Program Meet the minimal standards of the 1620 Users 

Group? Yes No_ 

Comment 

9. Were all necessary parts of the program received? Yes No_ 

Commen t 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 



Mr. Richard L. Pratt' Company 

Data Corporation 

7500 Old Xenia Pike Address 

Dayton, Ohio 45432 

^ User Group Code ( 



THIS REVIEW FORM IS PART OF THE 1620 USER GROUP ORGANIZATIONS PROGRAM 
REVIEW AND EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVITED 
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Simultaneous Equation Solution and Matrix Inversion With 
Complex Values 



Author: G. S. Haralampu 

New England Electric System 
441 Stuart St. 
Boston, Mass. 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 
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Revised - November 196k 



PROGRAM No* lit AND 701U - NEW ENGLAND ELECTRIC SYSTEM 
Ti'.Ie s Simultaneous Equation Solution and Katrix Inversion with Complex Values 
Written by t G. S # Haralampu, New England Electric S ye tea 

Scope s The program has been written for the solution of a maxima of fifteen 

equations with fifteen unknowns for the basic 1620, with the coefficients, 
unknowns, and constants being complex numbers* The matrix inversion method 
is used for the solution of the simultaneous equations* (Motet A 25 x 25 
matrix is allowed on the 7010 computer)* 

Machine Configuration s 20K 1620 with a 1622 card reader and punch, or 
80K 7010 computer 

Reference ; "Numerical Mathematical Analysis" by Scarborough* Fourth Edition, Page £U0* 

Speed : A 1* x h matrix requires 20 seconds for solution* A 16 x 16 matrix requires 12 
minutes for solution, on the 1620* The speed on the 7010 is approximately 
12 Urns s faster than the 1620* 

Analysis ; Assume the following h x k matrix is to be inverted : 



All the Z's are complex 



The following method was used: 

A. Place the first column of a unit matrix, equal to the siae of 
the main matrix, to the right of the main matrix. This is 
shown below t 

Main Matrix First Cola— of the Unit Matrix 



hi 


hi 


*13 


hk 


Z 21 


hi 


Z 23 




Si 


Si 


Z 33 




hi 


hi 


Z U3 


hk 



hi 


hi 


hj 


hk 


i 


hi 


hi 


Z 23 


hk 





Z 31 


Z 32 


Z 33 







hi 


hz 


Z U3 


hk 






B* By mathematical manipulation, obeying all the rules of matrix 
handling of rows and columns, obtain the first column of the unit 
matrix at the left of the main matrix* This will be: 
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First Column of Intermediate Values 

the Unit Matrix of the Main Matrix 



1 


*12 


*13 


m, 


*15 





*22 


*23 


*2U 


*25 





*32 


T 33 


*3U 


*35 





*l»2 


*U3 




*U5 



C. Shift the intermediate values of the matrix to the left by one column 
and place the second column of the unit matrix to the right. 

Main Matrix Second Column of the Unit Matrix 



*11 


*12 


*13 


111* 





*21 


*22 


Y23 


*2U 


1 


*31 


*32 


*33 


*3fc 







*U2 


*l*3 








D. Repeat "B B until we have: 

Second Column of Second Intermediate Values 
the Unit Matrix of the Main Matrix 






*12 


*13 




Y 15 


1 


122 


*23 


*2l» 


*2S 





*32 


*33 


*3U 


*35 





Tkl 


T U3 




*U5 



E. Repeat parts B and C until the matrix manipulation has been done for tht- 
four columns of the unit matrix. Then the values left in position 



*12 


*13 


111, 




T22 


*23 


Y2U 


T25 


*32 


*33 


Y3ii 


Y35 


*1*2 


TU3 







oompo3e the inverted matrix. 



Once the matrix is inverted, the unknowns can be found by matrix multi- 
plication. Assume the given equations to be solved are in the form of 



Equation 1 



where each bracketed symbol represents a matrix. Solving for the I's 
we have: 

[ij - [y] x [jj Equation 2 

where the Y matrix is the inverse of the Z matrix. 

If the equations to be solved are as shown below, 

z ll*l ♦ z 12*2 ♦ Z 13 I 3 + Z Ui T fc " E l 
z 21 x l + Z 22 J 2 + Z 23 I 3 + z 2i* I ii " E 2 
z 31*l * z 32*2 ♦ Z 33 J 3 + Z 3U I U * E 3 
^l 1 ! * Z U2 X 2 ♦ Z U3 I 3 ♦ Z UU J U " E U 
then Equation 1 expanded is as follows: 





h2 


Z U 


hk 










z 21 


2 22 


*23 


Z 2li 


X 






-\ 


hi 


Z 32 


z 33 


Z 3l. 




*3 








Z l»2 


Z U3 


2 W 











and Equation 2 is: 







"£l 


*12 


Yl3 


TU 


El 


12 




X21 


Y22 


*23 


^ X 


E 2 






^31 


Y32 


Y33 


Y 3U 


E3 






Y U1 


Y li2 


Y U3 


Y UU 


% 



Upon completing the matrix multiplation, the unknown I's are: 
l x - T U E! ♦ Y 12 E 2 ♦ Y 13 E 3 ♦ Y^ 

1 2 - Y 21 E! ♦ Y 22 E 2 ♦ Y23 E 3 + Y 2lA 

1 3 - Y 31 E 1 ♦ Y 32 E 2 ♦ Y 33 E 3 + Y 3 | i E ll 

X U ' Y a E l * Y l42 E 2 * Y U3 E 3 + Y Uh\ 
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USERS' GUIDE 
INPUT DATA FOR THE 1620 OR THE 7010 



First Card : Title of job with a minus sign in column 72. This will skip to a new 
page on the U07 far each separate solution. 

Second Card : (Right justified in the first five columns). The number of equations 
to be solved. (Fixed point) 

Subsequent Cards : (Floating point) The real part of the coefficient is entered in the 
first ten columns, and the imaginary part in the next ten colunns. 
For example, for a three equation solution, the Input 1st 



z ll 


z 12 


*L3 


El 


Z 21 


Z 22 


hi 


E 2 


Z 31 


z 32 


Z 33 


E 3 



and 12 cards are required* 

Last Card: (Fixed point; right justified in the first five colxmma) 

a) The value 1 If another case is to follow, and a 9999 if this was the 
last case* 



See 



iples in the Appendix. 



OUT POT 



LIMITATIONS 



1) Largest matrix handled is a 15 x 15 for the 1620, and 25 x 2$ for the 7010. 

2) No seroes allowed in the main diagonal* 

3) If a larger than a 15 x 15 on the 1620 (or a 25 x 25 on the 7010) Is attempted 

to be read in, the computer will not solve that case but will print message saying 
so, and then proceed to the next case* 
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FOR THE 7010 ONLY 

Place the following card ahead of the input data! 

(Columns 6-10) M0N$$ 

(Columns 16-28) EXEQ HARA,MJB 

A Matrix Inversion Subroutine has also been written for th- ?mn 
To make use of this subroutine do the following, written for the 7010. 

A. In your main program include the statements: 

1. DIMENSION R(25,26), X (25,26) 

2. COMMON R,X,K 

3. Read K or set K to equal size of matrix to be solved. 

i*. CALL MATTNV(R,X,K) 

B# l^° n nt °i thC ±tSpUt ' includc the following card before the »M0N$| 
EXEQ " card of your program: 

(Columns 16-26) CALL MATINV 
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OPERATING INSTRUCTIONS 



1) Clear Memory 

2) Load Object program in computer with the input for all the cases to be solved. 

3) When solution is completed for all the cases, the 1620 will PAUSE, and the 
7010 will EXIT. 



© 



INDEX TO THE APPENDIX 



Appendix I - Input and output from the 1620 

Appendix II - Fortran listing for the 1620 

Appendix III- Output from the 7010 

Appendix IV - Fortran listing for the 7010 

Appendix V - Fortran listing for the Matrix inversion subroutine for the 7010 

Appendix VI - Fortran program lietinf aolring aiamltaneoue equation*, using the 
■etrix inversion subroutine . This was written to test the MATIMV 
subroutine. 



mo 

Qo 

o 



Orn 



INPUT DATA 



APPENDIX I 



THIS IS A 4X4 MATRIX 



4 




0.071 


.363 


-8.023 


-4.031 


6. 146 


4.856 


_ _a«6*e 


6.323 


r.ooo 


0.0 


♦ 042 


•236. 


3.217 


7»501 


-3.31 1 


-7jl813_ 


• 665 


2.632 


1 .000 


0.0 


-t .452 


-5.067 


-9.569 


-9.686 


8*050 ~ 


4.390 


10. 167 


12.109 


1 .00 ~ 


0.0 


• 097 


• 564 


-8.101 


-4.564 


8.105 


4.718 


8.715 


7.022 


t .00 


0.0 



1 

THIS IS A 3X3 MATRIX 



3 

•0582 .168 



•0 -.00196. 

.0 -.000465 

♦058 ""' -.00502 
0»0 -.00196 
•05«i .168 
0.0 -.00196 



-.1102 .0000687 
0.0 -.000465 



0.0 -.00196 
•0582 *168 



-•1102 .0000687 
9999 



10 



57070 
57070 
57070 
57070 
57070 
57070 
57070 
57092 

57130 
57252 
57384 
57504 
57588 
57666 
57774 
57896 

57948 
57948 
57972 
58040 
58052 
58064 
58256 
58376 
58376 
68376 
58400 
58436 
58448 
58460 
58544 
58664 
58748 
58772 
58852 
58932 
59112 
59160 
09208 
09220 
59436 
59664 
59748 
59868 
59892 
59916 
59984 
T0020 
TO044 
T0052 
T0088 
Toil 2 
Tot 24 



APPENDIX II 

C PROGRAM NUMBER 14**N0V 1,1964 

C SIMULTANEOUS EQUATION SOLUTJON USING MATRIX INVERSION METHOD 
C COMPLEX NUMBERS ARE USED 

C WRITTEN BY G .S .HARALAMPU - NEW ENGLAND ELECTRIC SYSTEM 
C SIZE IS LIMITED TO A 15X15 MATRIX 

DIMENSION R(15,16),X(15,16),ER(15),EI(15),TR(15),TI(15) 

1000 FORMAT ( 1 5) 

1001 FORMAT ( 2F1 0.4,215) 

1002 FORMAT (/46H THE INVERTED MATRIX ELEMENTS ARE AS FOLLOWS) 

1003 FORMAT (// 46H THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW ) 

1004 FORMAT (48HMATR IX LARGER THAN 1 5X1 5 .NEXT CASE TO BE SOLVED.) 

1005 FORMAT (30H REAL I MAG ROW COL) 

1006 FORMAT (27^ REAL I MAG ELEMENT) 

1007 FORMAT (42H AN ELEMENT ON THE MAIN DIAGONAL IS ZERO) 

1008 FORMAT (49H ) 

1009 FORMAT(36X,35X,1H-) 
55 READ 1008 

C READ IN THE SIZE OF THE MATRIX 
READ 1000,K 
IF(K-16)103,101 ,101 
103 DO 50 J-1,K 
DO 30 J m 1, K 

30 READ 1001, R(I,J),X(I,J) 
50 READ 1001, ER(I),EI(I) 

C READ IN THE LAST CODE CARD. 1 FOR MORE CASES TO FOLLOW, AND 9999 
C FOR THE LAST CASE 

READ 1000.NCOOE 

KX « K + 1 

DO 33 I - 1 , * 

DO 48 L»1,K 

R(L, KX) - 0.0 
48 X(L, KX) - 0.0 

R(I,KX)«1.0 

IY - 1 

IF(R(l,1))34 f 31.34 

31 »F(X(t, 0)34,100,34 

34 T1«R<l f 1)*R(U1)+X{l f 1)*X(l ,1) 
TIR-R(M) 

TIX-X(M) 
DO 35 J-1.KX 

WR-(R( I , J)*T1R+X( I ,J)*T1X)/T1 
WU(X( I , J)*T1R-R( I , J)*T1X)/T1 
R(l, J) - WR 

35 X(l, J) - Wl 
IX -0 

IY - 2 

IF(I-K)37,38,37 

38 MX- 1-1 
MY-1 

GO TO 39 
37 MY m I + 1 
MX * K 

39 DO 40 L » MY, MX 



T0 160 
T0208 
T0256 
T0268 
T0484 
T0688 
T0844 
T1072 
T1096 
T1164 
T1244 
T1256 
T1268 
T1304 
T1448 
T1664 
T1700 
T1712 
T1760 
T1808 
T1820 
T2084 
T2336 
T2420 
T2576 
T2588 
T2600 
T2612 
T2624 
T2636 
T2648 
T2900 
T2912 
T2924 
T2936 
T3068 
T3 136 
T3148 
T3160 
T3172 
13180 
13192 
T3204 
13216 
T3288 
13360 
T3384 
T3392 

T9979 
T9949 
19909 
T9869 



40 



45 
44 



46 
33 



53 



49 



201 
202 
102 

100 

101 



2001 
2000 



T1 - R(L, 1) 
T2 « X(L, 1) 
DO 40 J«1 ,KX 

WR * R(l, J)*T1 - X(l, J)*T2 
Wl « R(l, J)*T2 + X(l, J)*T1 
R(L, J) - R(L, J) - WR 
X(L,J)-X(L,J)-WI 
IY - 3 

IF( 1-1)45,44,45 

IF((K-1)-IX)38,44,38 

DO 46 L - 1 , K 

DO 46 J - 1, K 

NU m J + 1 

R(L,J)«R(L,NU) 

X(L, J) - X(L, NU) 

CONTINUE 

DO 53 1-1, K 

TR(I )«0.0 

TKD-0.0 
DO 53 J-1,K 

T1-R(I,J)*ER(J)-X(I,J)*EI(J) 

T2«X(I ,J)*ER(J)+R(I,J)*EI(J) 

TR(I)-TR(I)+T1 

TI(I)-TI(I)+T2 

PUNCH 1009 

PUNCH 1008 

PUNCH 1002 

PUNCH 1005 

00 49 U1,K 

DO 49 J - 1, K 

PUNCH 1001, R(I,J),X(I,J),I,J 
PUNCH 1003- 
PUNCH 1006 
00 201 U!,K 

PUNCH 1001, TR(1 ),TI(I),I 

IF(NC0DE-1)102,55,102 

PAUSE 

PAUSE 

PRINT 1007 
GO TO 202 
PRINT 1004 
00 2000 U1, K 
00 2001 J-1,K 
READ 1001 ,W,Z 
READ 100t,W,Z 
REAO 1 000, N CODE 
GO TO 202 
END 



COS 
ATANF 
LOGF 
SIGNF 



S 



.9939 EXP 
.9899 SQRT 
9859 ABS 



T9999 SIN 

S9969 COSF 
9929 EXPF 
9889 SQRTF 
9849 ABSF 



[9989 SINF 
19959 ATAN 
[9919 LOG 
f9879 SIGN 
19839 R T7449 



IX 



IX +1 
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It 



T7U39 X 
1*599 



T4739 TR 
TUklS T001 
T4389 T005 
Tk3k9 1009 
TM09 BOO 
TU269 1 
T4229 *X 

WW B03* 
W09 BOt 

Jtofl S08 

T3989 NX 

If* I 2 ,. 
T3909 B046 

T3869 6201 

T3829 1001 

LOAD SUBROUTINES 



T5049 



T*A19 T002 
1*379 1006 
J4339 BOSS 
1*299 Sl03 
1*259 B03O 
1*219 BOOT _ 

f*179 B000000055 
4139 Kit 
TM)99 TtR 
f*059 5g2 
14019 BOO! 
13979 MY , 
13939 BOOl 
T3S99 JW 
T3859 6202 
13819 W 




T*889 El 
T**39 Tooo 
14399 lOO* 
1*359 T008 
1*319 BoiS 

8279 60S0 
239 MCOOE 

~*199 BOW 
*159 It 

179 6035 

^39 « ^ 
JW9 B038 
3959 >»*0 
3919 S>** 
3079 §0*9 
3839 1000 



T*7*9 



1 



REAL 
••1862 
-•0085 
-•0343 
.2154 

• 2298 

• 0297 
-•0082 
-.1673 

• 2177 
-•0080 
-•0058 
-.1386 

• 0)27 
•0197 
•0056 
.0175 



APPENTU X III 



THIS A 4X4 MATRIX 

THE INVERTED MATRIX E 
RO 



I MAG 
1.0951 

• 0151 

• 1262 
-1*2256 
-1.2040 

-•0935 
.0331 
1.1311 
-1.1513 
-.0049 

• 0339 
1.0672 
-•0725 
-.0779 
-•0046 

.0319 



EMEN 
CO 



OUTPUT 



ARE AS FOLLOWS 



THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW 
REAL I MAC CLEMENT 



-.0137 


• 0129 


I 


.0636 


-.1333 


2 


.0696 


-•0SS1 


3 


•0734 


-.1231 


4 



A p ? ,r rrax in 

OUTPUT 



THIS IS A 3X3 MATRIX 

THE INVERTED MATRIX ELEMENTS ARE AS FOLLOWS 



REAL 


I MAG 


ROW 


COL 


1.8*17 


-5.3149 


1 


1 


• 0385 


-.0488 


1 


2 


.0096 


-.0119 


1 


3 


.0385 


-.0488 


2 


1 


1.8*22 


-5.3153 


2 


2 


• 0385 


-.0488 


2 


3 


• 0096 


-.0119 


3 


1 


• 0385 


-.0488 


3 


2 


1.8417 


-5.3149 


3 


3 



THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW 
REAL I MAG ELEMENT 

.0748 -.3108 I 

-•2049 .5882 2 

-.2063 .5904 3 
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C PROGRAM NUMBER 14**MAY 1*1964 

C SIMULTANEOUS EQUATION SOLUTION USING MATRIX INVERSION METHOD 

C COMPLEX NUMBERS ARE USED 

C WRITTEN BY G. S.HARALAMPU - NEW ENGLAND ELECTRIC SYSTEM 

C SIZE IS LIMITED TO A 25X25 MATRIX 

DIMENSION R<25*26>*X<25*26> •ER(25)*EI (25) »TR<25)*TI <25) 

1000 FORMAT ( 15) 

1001 FORMATC2F10. 4*215) 

1002 FORMAT (/46H THE INVERTED MATRIX ELEMENTS ARE AS FOLLOWS) 

1003 FORMAT C///46H THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW ) 

1004 FORMAT (67H THIS MATRIX IS GREATER THAN 25X25* THE NEXT CASE WILL 
1 BE SOLVED. ) 

1005 FORMAT (30H REAL I MAG ROW COL) 

1006 FORMAT (27H REAL I MAG ELEMENT ) 

1007 FORMAT (42H AN ELEMENT ON THE MAIN DIAGONAL IS ZERO) 

1008 FORMAT < 1 X * 52H ) 
3000 FORMAT ( 1H1 ) 

55 READ (1*1 008 ) 
C READ IN THE SIZE OF THE MATRIX 

READ (1.1 000 ) K 
IF(K.GT.25)GO TO 101 
DO 50 IM«K 
DO 30 J • 1* K 
30 READ(1*I001) R(I*J)*XW*J> 
50 READ (1*1001) ER(I)«EI(I) 
C READ IN THE LAST CODE CARD. 1 FOR MORE CASES TO FOLLOW* ANO 9999 

C FOR THE LAST CASE 

READ < 1.1 000) NCOOE 
KX • K ♦ 1 
DO 33 I * 1 * K 
DO 48 L«l *K 
R (L * KX) ■ 0.0 
48 X(L* KX) • 0.0 
RM*KX)*l*0 
IY » 1 

IF(R( I* 1 ).NE.0.0)GO TO 34 
IF<X< 1*1 ).EQ.O.O >GO TO 100 

34 TI*R( I * i )*R< t*l)+X<I*l )*X( 1*1) 
TlR-R(I.l) 

T1X*X< 1*1) 
DO 35 J»l «KX 

WR«(R< I *J)#T1R+X( I • J)*T1X)/T1 
WI»CX< I • J)*T1R-R< UJ>*T1X>/Tl 
R < I • J > * WR 

35 XCI* <J> * WI 
IX «0 

IY » 2 

1F(I«NE.K)G0 TO 37 

38 NX- 1-1 
MY* I 

GO TO 39 
37 MY ■ I ♦ 1 
MX « K 

39 DO 40 L * MY* MX 
IX « IX +1 

Tl « R(L* 1 ) 
T2 ■ X(L* 1 > 



© 



00 40 JM*KX 

tfR m R(I, J)*T1 - X<I« J)*T2 
WI s R(I« J)*T2 + X<I» J)#T1 
R<L« J) n R(Li J) - WR 
40 X(L« J)«X<L« JI-WI 
IY « 3 

IF< I.EQ.l )G0 TO 44 

IF( (K-l )«NE*IX)GO TO 38 
44 DO 46 L • 1 • K 

DO 46 J ■ 1 • K 

NU « J ♦ 1 

R(L«J)*R(L«NU> 
46 X(Lt J> * X(L* NU) 
33 CONTINUE 

DO 53 I«ltK 

TR( I )«0«0 

TI ( I >«0«0 

DO 53 J«l#K 

Tl«R( I .J)«ER( J)-X(I « J)*EI t J) 

T2»X( I »J)«ER< J)+R< I »J)*EI < J) 

TR ( I )«TR( I )+Tl 
53 TI ( I )»TI ( I I+T2 

WRITE(3«3000) 

WRITE (3* 1006) 

WftITE(3«1002) 

WRITE (3tlO05) 

DO 49 I»1*K 

DO 49 J « I • K 
49 WRITE(3«100l ) P< I • J ) «X( I « J> • I 

WRITE <3« 1003) 

WRITE (3* 1006) 

DO 201 1-1 • * 

201 WRITEC3tl00l ) TR ( I ) tTI ( I ) ♦ I 

202 IF(NCO0E*EQ«1 )G0 TO 55 
102 CALL EXIT 

100 WRITE(3»1007) 
GO TO 202 

101 WRITE(3«l0O4) 
DO 2000 I»ltK 
DO 2001 J«ltK 

2001 READ (1*1001) W.Z 
2000 READ (1*1001) W.Z 

READ < 1 • 1 OOO ) NCOOE 

GO TO 202 

END 



© 



FORTRAN LISTING 1410-F0-970 
SUBROUTINE MAT INV 
DIMENSI0NA( 25,26) ,6(25,26) 
C0MH0NA,B*M 

01007 FORMAT ( 42H AN ELEMENT ON THE MAIN DIAGONAL IS ZERO ) 

KX»M*l 

00331*1, M 

D048L-1.M 

A(L,KX)=0.0 
00048 8UtKX>»0.0 

AU,KX)»L.O 

IV«i 

IF(AII, IUNE.0.0IG0T034 
IF(B< I, D-fcQ.O.OiGOTOlOO 

00034 Tl»Al I,1)»A( I,l)*B<I,l>»B< I, l> 
TiR*A(I f l» 

TlX»BU»i> 
D03!>J»1#KX 

WR»(A< I ,J)»UR*B( I, J)»T1X»/Tl 
Wi*(BU t J>*TlR-All,J>»TlXI/M 
AU,J) = HR 

00035 B(I,J)=WI 
IX*0 
IY»2 

IF( l.NE.M)G0T037 

00038 NX«I-i 
MY*1 
G0T039 

00037 MY»I*1 
MX*M 

00039 D040L*MY,MX 
IX*IX*i 
Tl*A(L,i> 
T2*B(Lt I) 
D040J«i,KX 

WR-AU, J»»T1-BU. JI»T2 
WI*A«I*J)»T2+BIItJ)*Tl 
A(L*JI*A(Lt Jl-WR 

00040 B(L, J)»B(L, J»*WI 
IY«3 

IF ( I «EQ« 1 ) GOT 044 
IFUM-1 >.NE,IX»G0T038 
00044 D046L«l,M 
D046J-I.M 
NU»J*i 

AILtJ)«A(LtNUI 
00046 8(L»J)-B(L,NU) 
00033 CONTINUE 

RETURN 
00100 HRITEI3, 1007) 

STOP 

END 



APPENDIX VI 



FORTRAN LISTING 1410-FO-970 
OINENS10NRt25,26J,XC25.26).ERC25l,fcIC25>»TR(25), TIC 25) 
COHHONRfXtK 

01000 FORMAT! 151 

01001 F0RMATC2F10.4*2I5> 

01002 FORMAT! M6H THE INVERTED MATRIX ELEMENTS ARE AS FOLLOWS* 

01003 F0RMAT!///46H THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW > 

01004 FORMAT I 6 7H THIS MATRIX IS GREATER THAN 25X25. THE NEXT CASE MILL B 

01005 FORMAT (30H REAL I MAC ROM COL I 

01006 FORMAT! 2 7« REAL I MAG ELEMENT I 

01007 FORMAT 1 4 2H AN ELEMENT ON THE MAIN 01 AGONAL IS ZERO I 

01008 FORMAT t IX, 52H > 
03000 FORMAT ClHi) 

00055 READ! 1» 1006) 

C REAO IN THE SIZE OF THE MATRIX 

READC1.10001K 

IFU.GT.251G0T0101 

DOSOI-l.K 

D030J«ltK 
00030 REAOC It 1001 )R< It J).X(I.J» 
00050 REAOC It 1001 )ERC I ) t El C I ) 

C REAO IN THE LAST CODE CARD. I FOR MORE CASES TO FOLLOW, ANO 9999 

C FOR THE LAST CASE 

REAOC It 1000 )N€ ODE 

CALLMATIMV 

00531-1. K 

TRCI>*0.0 

Til 11*0.0 

0053 J«UK 

Tl«Rf t, J MERC J>-XCItJMEIC Jl 
T2»XC I,J)*ERC J)+RCItJ>*EICJ> 
TRCIt«TRtI»+Tl 
00053 TIC I )»TI I I »^T2 
WRITEC3.3000) 
WRITE<3tl008> 
WRITEC3,1002) 
WRITfcC3tl005t 
00*91-1. K 
0049J«1,K 

00049 MRITEC3,100l»R(l,J).X(I,JI,I,J 
MAITEC3. 10031 
WR1TEC3.1006) 
002011- I, H 

00201 WRITEf3tl001)rRCIItTICn 9 1 

00202 IFCNCQ0E.EQ.1IG0T055 
00102 CALLEXIT 

00100 WRITEC3.1007) 
G0T0202 

00101 WRITEI3,1004) 
G0T0102 

END 



OUTPUT - CASK 1 



THIS IS A 4X4 MATRIX 

THE INVERTED MATRIX ELEMENTS AWE AS FOLLOWS 



EAL 


I MAG 


ROW 


COL 


- . 1 862 


1 .0951 


1 




-.0085 


.0151 


1 


2 


-•0343 


• 1282 


1 


3 


• 2154 


-1.2256 


1 


4 


• 2296 


-1 .2040 


P. 


1 


• D297 


-•0935 


2 


2 


-.O062 


.0331 


2 


3 


-.1673 


I .131 1 


2 


4 


• 2177 


-1 • 151 3 


3 


1 


-.0080 


-.0049 


3 


2 


-•0056 


.0339 


3 


3 


-•1386 


1 »0672 


3 


4 


• 0327 


-•0725 


4 


I 


• 0197 


-•0779 


4 


2 


.0054 


-•0046 


4 


3 


• 01 75 


.0319 


4 


4 



THE SIMULTANEOUS EQUATION SOLUTIONS FOLLOW 
REAL I MAG ELEMENT 

-•0137 .0129 1 

•0858 -.1333 2 

.0654 -.0551 3 

•0754 -.1231 4 
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OUTPUT - CASE 2 



THIS IS A 3X3 MATRIX 

THE INVERTED .MATRIX ELEMENTS ARt AS FOLLOWS 



REAL 


I MAG 


ROW 


COL 


1 .84 1 7 


-5.3149 




1 


• 0385 


-•0488 




' ?. 


• 0096 


-•0119 




3 


• 0385 


••0488 


2 


I 


t .8432 


-5.3153 


2 


? 


• 0385 


-.0488 


2 


3 


• 0096 


-.0119 


3 


I 


.0385 


-.0488 


3 


2 


1 .841 7 


-5.3149 


3 


3 



THE SIMULTANEOUS EQUATION SOLUTIONS EOLLO* 
REAL I MAG ELEMENT 

.0748 -.3108 1 

-.2049 .5882 2 

-•2063 .5904 3 



COMPUTER 
TECHNOLOGY 



1 026 2034 8 



